Determination of critical current density from arbitrary flux relaxation process 
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The current-carrying ability of a type-II superconductor is generally represented by its critical 
current density. This can be determined by measuring a flux relaxation process starting with a 
testing current density that is greater than or equal to the critical value. Here we show that a flux 
relaxation process starting with an intermediate current density can be converted into a process 
starting with the critical current density by introducing a virtual time interval. Therefore, one may 
calculate the critical current density from the flux relaxation process starting with a current density 
below the critical value. The exact solutions of the time dependence of current density in the flux 
relaxation process were also discussed. 

PACS numbers: 74.25.Wx, 74.25.Uv 
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I. INTRODUCTION 

Critical current density is the maximum current den- 
sity that can be carried by a type-II superconductor. Its 
value is mainly determined by the vortex pinning poten- 
tial. One can determine the pinning potential from a 
flux relaxation process, which is the phenomenon that 
quantized vortices jump between adjacent pinning cen- 
ters spontaneously due to various reasonsi"^. The flux 
relaxation phenomenon can be described by the Arrhe- 
nius equation in case of thermal activation. This suggests 
that a flux relaxation process has a strong dependence 
on temperature and vortex activation energyj^ The vor- 
tex activation energy is a decreasing function of current 
density because it can be reduced by the Lorentz force of 
the current density. A number of current dependent vor- 
tex activation energies^"— were proposed on the basis of 
different physical considerations. Using these vortex ac- 
tivation energies, one can obtain the corresponding time 
evolution equations of the current density.^ 

However, the time evolution equations of the current 
density give the critical current density at zero time. This 
means that, to find out the critical current density of a 
type-II superconductor, one has to apply a current den- 
sity that is greater than or equal to the critical value. But 
at lower temperatures, the critical current density may be 
largei^ and the measurements are technically difficult. In 
fact, many measurements were carried out with an initial 
current density well below the critical valucii^""— Thus, it 
is desirable to find a general formula which can be used to 
calculate the critical current density by measuring a fiux 
relaxation process with an initial current density below 
the critical value. 

For the above reasons, one should consider the vortex 
pinning potential, which plays a dominant role in deter- 
mining the critical current density of a type-II supercon- 
ductor. It is known that, at a constant temperature, the 
vortex pinning potential does not change in a fiux relax- 
ation process although the vortex activation energy does. 
This indicates that the effects of the vortex pinning po- 
tential are embodied in each vortex jump. Therefore, a 
vortex system may carry the information about the crit- 



ical current density even with a testing current density 
well below the critical value. An early study has shown 
that the maximum internal field in a vortex penetration 
or flux relaxation process can be calculated by introduc- 
ing some time parameters. In case of current density, it 
should be possible to construct a theoretical model with 
which one can predict the critical current density by mea- 
suring a flux relaxation process with an arbitrary initial 
current density. 

In this work, we first studied the time dependence of 
the vortex activation energies in two thought fiux relax- 
ation processes. The mathematical expression and physi- 
cal meaning of a virtual time interval are discussed. Next, 
we used the inverse-power activation energy and loga- 
rithmic activation energy to construct the time evolution 
equations of the current density. Finally, we discuss the 
exact solutions of the time dependence of current density 
with linear activation energy and logarithmic activation 
energy. Let us now derive the generalized expression for 
the time dependence of the vortex activation energy. 



II. TIME DEPENDENCE OF ACTIVATION 
ENERGY 

In the flux relaxation process, the current density j 
can reduce the vortex activation energy Ua by its Lorentz 
force and therefore increase the vortex hopping rate. This 
indicates that Ua is a decreasing function of j. Because j 
is a decreasing function of the time t, Ua is an increasing 
function of t. 

The rate of change of the current density, dj/dt, is pro- 
portional to the vortex hopping rate. Furthermore, the 
current density in a flux relaxation process is a decreas- 
ing function of the time, i.e., dj/dt < 0. According to 
the Arrhenius equation, therefore, we have the following 
equationi^""— 



dt 



-T- = ~Ce 



-Ua/kT 



(1) 



where C is a positive proportional constant. 

The exact solution of Eq.(IT]) is currently unavailable. 
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but we can find an approximate solution to it. Consider 
a superconductor in which an arbitrary current density 
starts to decay at zero time {t = 0). The initial vortex 
activation energy is Ui. As the time increases to t, the 
vortex activation energy increases to Ua- Rewrite Eq.([T]) 
and integrate it on both sides 



e^^/^^dUa = - / C^dt. 



With logarithmic accuracy^, we obtain the following 
equation 



(2) 



where r = —kT/[C[dUa/dj)\ is a short time scale 
parameter—. Here we included a negative sign to ensure 
that r is a positive number because in a flux relaxation 
process the activation energy Ua is a decreasing function 
of the current density j, i.e., dUa/dj < 0. 

Eq. ^ is the general time evolution equation of the vor- 
tex activation energy in a flux relaxation process start- 
ing with an arbitrary current density. If the flux relax- 
ation process starts with the critical current density jc, 
then the initial vortex activation energy Ui is zero, i.e., 



U^ 



0. Therefore, Eq.® reduces to 



Uacit) = kTln ( 1 + - 



(3) 



Eq.Q is the well-known time dependent vortex acti- 
vation energy used in earlier theories^^, from which the 
current decay models were derived. 

Because our purpose is to calculate the critical current 
density jc from a flux relaxation process with an initial 
current density ji below the critical value, we need to 
consider the possibility of converting Eq. ^ into the form 
of Eq.Q. To do this, let us first study the following two 
thought flux relaxation processes; 

Case A: At temperature T, a vortex system in a su- 
perconductor starts to decay with an initial current den- 
sity ji below the critical value where the vortex activa- 
tion energy is Ui. The time dependence of the vortex 
activation energy is then governed by Eq.Q, that is, 
Ua{t)^kTlnie^^/'''^ + t/T). 

Case B: At temperature T, consider the same super- 
conductor used in Case A. But a vortex system starts 
to decay with the critical current density jc, where the 
initial vortex activation energy is 0. After a time in- 
terval ti, the vortex system and current density reduces 
to exactly the same as that in Case A, where the vor- 
tex activation energy is Ui. If we now restart count- 
ing the time t, then the flux relaxation process should 
be governed by Eq.dS]) with an argument ti + t, that is, 
Uacik + t) = kTln[l + {U + t)/T]. 

Since Case A and Case B describe the same physical 



process at time t, they must obey the same physical law. 
Therefore, the vortex activation energies should be equal 
in quantity for all time t, that is. 



Uait) = UaciU+t). 



(4) 



Substituting the exact expressions of Ua{t) and 
UaciU + t) into Eq.®, we have 



(5) 



Physically, we can consider the decay process during 
the time interval ti in Case B as a "virtual process" in 
Case A. This is equal to doing a time transformation 
t' ^ ti +t in Eq.®. Therefore, Eq.© can be rewritten 
as 



Uait) = kTln 1 + 



t,+t 



(6) 



where ti is defined by Eq.([5]). Substituting Eq.® into 
Eq.®, it goes back to Eq.®. This means that the phys- 
ical considerations on ti is consistent in mathematics. 

It is now clear that, by introducing the virtual time 
interval ti, the vortex activation energy of a flux relax- 
ation process starting with a lower current density ji can 
be converted into that of a process starting with the crit- 
ical current density jc. The physical meaning of ti is: 
the virtual time interval during which the vortex activa- 
tion energy Ua{t) increases from to the initial value Ui. 
This will be further discussed later in the study of the 
time dependence of current density. 



III. TIME DEPENDENCE OF CURRENT 
DENSITY 



The time dependence of current density j{t) can be 
obtained by combing the time dependence of the vor- 
tex activation energy Ua{t) (see Eq.®) with the current 
dependence of the vortex activation energy Ua{j). Var- 
ious Ua{j) were proposed and each has particular phys- 
ical meaning. ^^^^ Here we used the inverse-power acti- 
vation energy^ and logarithmic activation energy^i^ to 
construct the j{t) respectively. 



A. Inverse-power activation energy 

In case of random pinning, collective creep theorySil 
gives a good description to the vortex behavior. This the- 
ory gives a vortex activation energy that has an inverse- 
power dependence on current density, that is, 



UaU) 



Uo 



1 



(7) 



where /i is a critical exponent related the vortex structure 
in the superconductor. The number 1 on the right side 
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is introduced to ensure that Ua{jc) = 0. 

Substituting Eq.® into Eq.©, we have the time de- 
pendence of the current density 



j{t) = jc 



1 + -Tj—ln I 1 + 



(8) 



Putting t = in Eq.®, we obtain the initial current 
density 



Ji — Jc 



1 



In 1 H 

\ T 



1 -i/a* 



(9) 



Fitting Eq.® to experimental data, we can obtain the 
critical current density jc (see Fig. 1(a)). This means that 
we can calculate jc by measuring a flux relaxation process 
starting with an arbitrary current density ji. This oper- 
ation is possible because jc is mainly determined by the 
vortex pinning potential, which affects each vortex jump 
(see Eq.([T])). It does not matter what current density the 
flux relaxation process starts with. 

Eq.® also shows that, by introducing the virtual time 
interval t^, a flux relaxation process starting with a cur- 
rent density ji can be converted into a process starting 
with the critical current density jc- This can be analyzed 
in a way analogue to that used in Case A and Case B. 
The process during U is a virtual process, introduced for 
assisting physical understanding and easy calculation. 

In Eq.([S]), we expressed the virtual time interval ti in 
terms of the initial activation energy Ui. But in practice 
the measurable physical quantity is the initial current 
density ji. Therefore, it should be more convenient to 
express ti in terms of j^. Inverting Eq.®, we have 



exp 



Up 



(10) 



Putting ji — jc in Eg. ((TO)) , we have ti = 0. Therefore, 
the physical meaning of the virtual time interval ti can 
also be explained using the current density: the time in- 
terval during which the current density reduces from the 
critical value jc to the initial value ji . 

Fig. 1(a) shows the fitting of Eq.® to the experimen- 
tal data from a Bi2Sr2CaCu20s+x single crystal. The 
parameter /i = 0.12 indicates that the vortices in the 
sample are approximately "single vortex" ^ 



B. Logarithmic activation energy 

In high-Tc superconductors, sometimes it is more accu- 
rate to use the following logarithmic activation energj«2iS 



Ua{j) = Uoln ( ^ 



(11) 



Substituting Eq. ([5]) into Eq. dTTI) , we have the time de- 
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FIG. 1: (Color online) Time dependence of current density. 
The scattering points are the persistent current density in- 
duced in a Bi2Sr2CaCu20s+x single crystal at 26 K. The 
solid black lines are the theoretical fits, (a) Fitted with the 
inverse-power model (Eq.©). jc = (1.22 ±0.01) x 10^ A/m'^, 
At = (0.12 ± 0.01), Uo = (79.8 ± 0.1)fc, ti = (0.81 ± 0.01) s, 
r = (1.20 ± 0.01) s. (b) Fitted with the logarithmic model 
(Eq.dlll)). jc = (1.20 ±0.00) X 10**A/m^ Uo = (92.2 ± 0.1)fc, 
ti = (0.79 ± 0.01) s, r = (1.03 ± 0.01) s. 



pendence of the current density 



ti+t 



]{t)^Jc 1± 



'kT/Uo 



(12) 



Putting t = in Ea.(IT2)). we obtain the initial current 
density. 



ji = ic 1 1 ± - 



-kT/Uo 



(13) 



Inverting Ea.(|13p. we obtain the virtual time interval 

Uo/kT 



ti = T 



1 



(14) 



Fig. 1(b) shows the fitting of Ea. ([T^ to the experimen- 
tal data from a Bi2Sr2CaCu20g,+x single crystal. 



IV. EXACT SOLUTIONS OF THE TIME 
DEPENDENCE OF CURRENT DENSITY 

We have used Eq.® in the derivation of Eq.® and 
Eq. ([T2|) . Eq.® provides a universal method for calculat- 
ing the time dependence of current density with various 
vortex activation energies. But the solutions obtained in 
this way are approximate solutions. Here we show that 
one can find an exact solution of Eq.([T|) using the linear 
activation energy Ua{j) = Uo{l — j/jc) and logarithmic 
activation energy Uaij) = Uoln{jc/j). 



A. Linear activation energy 

Substituting the linear activation energy Ua{j) 
Uo{l — j/jc) into Eq.dl]), we have 



dt 



= _(^g-/3(l-j7jc) 



(15) 



where /S = Uo/kT. 

The solution of Eq.din]) is 



j{t) = Jc 







1 - ii„ (l + 









(16) 



where t = jc/ {f3C). 

Putting t ~ in Eg. (1161) . we obtain the initial current 
density, 



=Jc 



(17) 



Inverting Eg. (1171) . we obtain the virtual time interval 

B. Logarithmic activation energy 

Substituting the logarithmic activation energy Ua{j) = 
Uoln{jc/j) into Eq.([T]), we have 



dt \]c 



The solution of Eq.(|T9l) is 



j{t)=]c 1 + 



U + t 



1/(1-/3) 



(19) 



(20) 



where t = jc/[{fi - 1)C]. 

Putting i = in Eg. (l20l) . we obtain the initial current 
density, 



j^ = Jc 1 1 + - 



1/(1-/9) 



(21) 



Inverting Eg. (j2ip . we obtain the virtual time interval 

(1-/3) 



U = T 



1 



(22) 



V. DISCUSSION 



Eq.®, Eq.(lT2]), Eq.([T6l), and Eq.(|20l) are the time evo- 
lution equations of current density with generalized ini- 
tial conditions. Using these equations, one can calculate 
the critical current density jc from a flux relaxation pro- 
cess with a lower initial current density ji . 

The inverse-power activation energy Eq.([7]) is obtained 
from the collective creep thcory^'^, which has profound 
physical meaning. This model also gives the information 
about the vortex structure by the exponential parame- 
ter n. Therefore, one may choose this model if he/she 
is interested in finding out the vortex structure in the 
superconductor. 

Sometimes, it is found that the logarithmic activation 
energy Eg. ljlip gives a better description to experimen- 
tal data^^i^ The corresponding time evolution equation 
of the current density Ea. (|12p has four free parameters. 
Therefore, this model is more convenient for analyzing 
experimental data. One may choose this model if he/she 
is only interested in finding out the current dependence 
of the vortex activation energy U{j). 

At higher temperatures, however, the flux relaxation 
phenomena usually exhibit strongly non-logarithmic be- 
havior. The experimental data usually cannot be well 
fitted by the theoretical models. In this case, one may 
consider using the infinite series model^iii^, which gives 
the information about the elastic and non-elastic defor- 
mation of the vortices. This generalized model includes 
the linear model (see Eg.lfTSI) as a special case. Using 
this theory, one can calculate the activation energy with 
high accuracy. But this theory also includes many free 
parameters and one has to choose the higher order terms 
on his purpose. 

Finally, it should be emphasized that in Ref . [Ill we pro- 
posed a generalized infinite series expression for the cur- 
rent dependence of the vortex activation energy, which is 
related to the difference of the superconductivity in var- 
ious superconductors. This enables us to analyze the ex- 
perimental data on the current decay behavior and then 
calculate the activation energy of a vortex system with- 
out subjecting to any special conditions. In the present 
work, however, we introduced the virtual time interval 
ti for the time dependence of a flux relaxation process. 
This is equal to constructing a generalized initial condi- 
tion for the time dependence of a flux relaxation process. 
It is related to the "time" . 
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VI. CONCLUSION 

The information about the critical current density of a 
type-II superconductor is carried by the vortex system 
throughout a flux relaxation process. By introducing 
a virtual time interval, a flux relaxation process start- 



ing with a lower current density can be converted into 
a process starting with the critical current density. This 
enables us to compute the critical current density from 
a flux relaxation process with a lower current measure- 
ment. It does not matter what initial current density the 
flux relaxation process starts with. 
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